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The antiferromagnetic Heisenberg model is studied on a two-dimensional bipartite quasiperiodic 
lattice. The distribution of local staggered magnetic moments is determined on finite square ap- 
proximants with up to f393 sites, using the Stochastic Series Expansion Quantum Monte Carlo 
method. A non-trivial inhomogeneous ground state is found. For a given local coordination num- 
ber, the values of the magnetic moments are spread out, reflecting the fact that no two sites in a 
quasicrystal are identical. A hierarchical structure in the values of the moments is observed which 
arises from the self-similarity of the quasiperiodic lattice. Furthermore, the computed spin structure 
factor shows antiferromagnetic modulations that can be measured in neutron scattering and nuclear 
magnetic resonance experiments. This generic model is a first step towards understanding magnetic 
quasicrystals such as the recently discovered Zn-Mg-Ho icosahedral structure. 



PACS numbers: 

Quantum magnetic phases of low-dimensional antifer- 
romagnetic (AF) Heisenberg systems are known to show 
various degrees of disorder caused by zero-point fluctua- 
tions. For example, isolated spin-1/2 chains have quasi- 
long-ranged antiferromagnetic order |Q , whereas two-leg 
ladders are short-range-ordered and true long-range 
order is found in the two-dimensional square lattice 
While these states are homogeneous due to the trans- 
lational invariance of the underlying lattice, it is inter- 
esting to explore how non-periodic environments, such 
as provided in quasicrystal structures, affect the mag- 
netic properties of quantum magnets. Naturally, it can 
be expected that instead of a uniform order parameter, 
such as the staggered magnetization in a Heisenberg an- 
tiferromagnet, there will be instead a distribution of lo- 
cal order parameters. The magnitude of each local mo- 
ment may, for example, strongly depend on its environ- 
ment. In a local viewpoint, the magnetic moments on 
sites with a greater number of neighbors are expected to 
be suppressed due to increased local spin fluctuations. 
|0| This picture will be modified by the quasicrystalline 
equivalent of non-local spin- wave excitations, leading to a 
non-trivial distribution of the order parameters. Recent 
inelastic neutron scattering experiments on the Zn-Mg- 
Ho isocahedral quasicrystal || have revealed an antifer- 
romagnetic superstructure, which fits very well with the 
antiferromagnetic Heisenberg model. As opposed to mag- 
netic quasicrystals with itinerant charge carriers 0, the 
electrons appear well localized in this material, and dis- 
play an antiferromagnetic modulation with a large wave 
vector at temperatures below 6 K, similar to the pattern 
of the generic quasicrystal structure reported here. 

In this work we use the recently developed Stochas- 
tic Series Expansion Quantum Monte Carlo (SSE-QMC) 
algorithm [Q to analyze the magnetic ground state prop- 
erties of the nearest-neighbor AF spin-1/2 Heisenberg 
model, 



H = jJ2^i-S j , (1) 

on a two-dimensional quasiperiodic tiling. This bipartite 
structure, shown in Fig. 1, is called "octagonal tiling" 
due to its overall 8-fold symmetry. Sites in this tiling 
have coordination numbers z ranging from 3 to 8. This 
generic lattice structure was chosen to ensure that the 
magnetic ground state is unfrustrated due to the bipar- 
tite property, i.e. nearest-neighbor sites belong to two 
distinct subtilings. j^] In the numerical study, approx- 
imants with N = 41,239, and 1393 sites || are con- 
sidered, and toroidal boundary conditions are applied. 
ifiof The temperature is chosen low enough to obtain the 
ground state properties of these finite systems. 




FIG. 1. Approximants of the octagonal tiling, (a) Inflation 
from 239 (thin lines) to 41 (thick lines) sites, (b) "Heisenberg 
stars" are the local building blocks. 

An important property of quasicrystals, in the absence 
of invariance under translations, is their self-similarity 
under inflation transformations. These are site decima- 
tion operations that increase the length scale but glob- 
ally preserve the quasiperiodic structure. An example of 
inflation is shown in Fig. 1(a) where original tiles and 
inflated tiles are superimposed. Since periodic approxi- 
mants are finite pieces of the octagonal tiling, they are 
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not invariant under such a transformation, but instead 
transform into smaller approximants. In this process, 
the sites which remain after an inflation operation can 
change their connectivity. A site of coordination number 
z will be transformed into a site of new coordination num- 
ber z' after inflation, with z' < z. We are interested in 
how this self-similarity is reflected in the inhomogeneous 
magnetic ground state of the Heisenberg Hamiltonian. 

In the SSE-QMC simulations the local value of the 
staggered magnetic moment is determined at each lattice 
site i. It is defined by 
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where the sum extends over all lattice sites j of the 
approximant. The inset of Fig. 2 shows a finite-size 
extrapolation of the spatially averaged staggered mag- 
netic moment, m = ^(3/A 2 ) )> wnerc 
the sum extends over all pairs of sites i and j. In the 
thermodynamic limit, this observable approaches m = 
0.337 ± 0.002, indicating that the system has AF long- 
range order. Note that the average moment is larger than 
that of the square lattice, m > m s = 0.3071 ±0.0003 [jfl), 
even though the average connectivity for the octagonal 
tiling is exactly z = 4, same as for the square lattice. 
This suggests that quantum fluctuations reducing the or- 
der parameter are suppressed due to the inhomogeneous 
connectivity in the octagonal tiling. 
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FIG. 2. Dependence of the staggered magnetic moments on 
the coordination number z. Inset: finite-size extrapolation of 
the average staggered magnetic moment. 

The measured dependence of m(i) on the local coor- 
dination number z is plotted in the main part of Fig. 2 
for the N=1393 approximant. A wide spread of these 
moments is observed, particularly for small values of z. 
There is a tendency for sites with more neighbors to have 



smaller staggered moments because their spin is sup- 
pressed by forming a larger number of local bonds. Q] 
However, there are exceptions to this rule, indicating the 
presence of longer-range correlations. A further striking 
feature observed in Fig. 2 are splittings in the distribu- 
tion function for the staggered moments. This is most 
evident at z = 5, but it also occurs on smaller scales at 
other coordination numbers, e.g. at z = 8. 

Local environments clearly play an important part. To 
sort out the roles of the local versus the long-range corre- 
lations, we next focus on the local approximation to the 
estimator for the staggered moment, 



mioc(i) 
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z + 1 
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where the sum is truncated to the sites in the immediate 
vicinity of i. This quantity is a local approximation of 
m(i) defined in Eq. (2). It reflects the average bond 
strength (1/z) YTj=i ^« ' Sj of site i with its z neighbors. 
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FIG. 3. Dependence of the local staggered magnetic mo- 
ments on the coordination number z. The points on the 
dashed line are exact values obtained for Heisenberg stars. 
Inset: spatial variation of the local approximation to the local 
staggered moments where the radii of the circles correspond 
to the size of local staggered moments. 

Fig. 3 shows a clear monotonic decay in mi oc (i) with 
increasing number of neighbors z. This effect can be 
understood by considering the local constituents of the 
cluster, shown in Fig. 1(b), which we call Heisenberg 
stars (HS). These are described by the nearest- neighbor 
Hamiltonian H HS (z) = J^2* =1 So ■ Sj, where the index 
denotes the central spin of the cluster and the AF ex- 
change integral is set to J = 1. The ground state energy 
of these stars can be calculated exactly, and is given by 
Eq S (z) — — (2 + z)/4. The estimator for the local mag- 
netization can be obtained directly from Eq. (3), 
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implying a z _1 / 2 -like decay, consistent with the trend of 
the QMC data observed in Fig. 3. The ground state 
energy as well as the magnetic moment of the isolated 
Heisenberg clusters are naturally larger than the equiva- 
lent quantities in the quasicrystal where the presence of 
additional bonds lowers the energy per site and further 
screens the local magnetization. 

For all coordination numbers, there is some degree of 
spread in the distribution of mi oc {z). This phenomenon 
can be understood by considering non-local corrections to 
the magnetic moments of the isolated Heisenberg stars. 
Next-nearest-neighbor corrections to the local moments 
for the six types of sites (z — 3, 8) can be obtained in a 
self-consistent manner by re-diagonalizing the Heisenberg 
stars, using local spin operators that are renormalized ac- 
cording to their specific environments, i.e. by replacing 
(Sf) by mf^c (z) with the appropriate coordination num- 
ber to determine the renormalized matrix elements. This 
procedure naturally leads to larger spreads for sites with 
small z, since in the octagonal tiling they have a wider 
variety of local environments than the sites with large z, 
which is observed in Fig. 3. 

Furthermore, one notices a prominent bimodal split- 
ting effect at sites with z = 5. In the octagonal tiling, 
these sites occur in pairs and have two distinct local en- 
vironments. |H| The first type of z = 5 clusters has two 
z = 4 neighbors, two z = 3 neighbors, and one z = 5 
neighbor. The other type has four z = 4 neighbors and 
one z = 5 neighbor. Analogous to the discussion in the 
above paragraph, the next-nearest-neighbor correction to 
the staggered magnetic moment for the two types can be 
determined by diagonalization of the corresponding 6x6 
matrices in the Sf ot — (N — 2)/2 subspace of these two 
N=6-site Heisenberg stars, using the magnetization val- 
ues (z) of the Heisenberg stars in Eq. 4. Jl3| Follow- 
ing this procedure, one finds two renormalized staggered 
moments with a splitting Ato/ oc (5) = 0.007, which is 
in very good agreement with the numerically observed 
splitting for the z — 5 sites in Fig. 3. 

Let us now turn to the hierarchical structure of the 
magnetic ground state that can be observed in the split- 
tings of the staggered moments of the z = 8 sites. These 
sites have eight-fold symmetry out to their third neigh- 
bor shell. Since their local environments are identical up 
to a larger distance compared to other sites, the spread 
in values of mi oc (i) is relatively small (Fig. 3). One 
can however divide these sites into four groups accord- 
ing to their transformation under inflation, i.e. their new 
local coordination numbers z' . The m(i) values conse- 
quently fall into four discrete groups [z = 5, 6, 7, and 8), 
as shown in Fig. 4. This hierarchical fine structure is 
analogous to the discussion of the z = 5 sites, but occurs 
on a smaller scale. Interestingly, the most symmetric 



group with z' = 8 experiences an additional hyperfine 
splitting due to a four-fold subgrouping with z" = 5, 6, 7, 
and 8 under further inflation. In the limit of infinite size, 
this leads to a multifractal distribution of m(i) for this 
class of sites. 
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FIG. 4. Hierarchy of the staggered moment of the z — 8 
sites in the 1393 site tiling, grouped according to the value 
of z' under the inflation transformation. Numbers on top of 
the symbols give the value of z" for the z — 8 sites under a 
further inflation transformation. 

To summarize this discussion, the different local envi- 
ronments in the quasiperiodic structures thus lead, firstly, 
to a systematic decrease of m/ oc (i) with increasing z. The 
differences in next-nearest neighbor shells for a given z 
gives rise to a spread of m/ oc (i) and to discrete distribu- 
tions which are self-similar on smaller and smaller scales, 
as seen in the 8-fold site example. These generic features 
should be observable by high-precision nuclear magnetic 
resonance measurements. Comparing the data in Figs. 2 
and 3 one further observes that long-range correlations 
tend to suppress the differences in m/ oc (i) between the 
sites, leading to a non-monotonous dependence of m(i) 
on z. The spread in values of m(i) is overall smaller 
than for mi oc (i). However, the splittings at z = 5 and 
z = 8 persist in m(i), and should thus be experimentally 
observable. 

We conclude the discussion of the magnetic ground 
state on the octagonal tiling by presenting the magnetic 
structure factor, which is of relevance to inelastic neu- 
tron scattering experiments. It is obtained by a Fourier 
transformation of the real space correlation function of 
the largest available approximant, 

^(k) = i^e ik -(- r ')(5 r z ^) ! (5) 

r,r' 

using the QMC data for the octagonal tiling. In Fig. 
5 the diffraction pattern of this quasicrystal lattice is 
compared with the antiferromagnetically induced super- 
structure. The eight-fold pattern in Fig. 5(a) reflects the 
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eight-fold symmetry of the quasicrystal. The introduc- 
tion of AF correlations leads to a splitting of Bragg peaks 
into the largest available wave vectors at the Brillouin 
zone boundary and their symmetry points, as it is seen 
in Fig. 5(b). For example, the central Bragg peak splits 
into Q w (±7r, ±7r/2) and (±7r/2,±7r), consistent with 
the modulation expected for the octagonal tiling. This 
AF superstructure is an extension to the two-dimensional 
case of the AF ID structure discussed in pj . 
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FIG. 5. (a) Simple Bragg scattering, and (b) antiferromag- 
netic superstructure in the octagonal tiling, extracted from 
Quantum Monte Carlo data on the 1393 site approximant. 

In conclusion, we have examined the AF spin- 1/2 
Heisenberg model on the bipartite octagonal tiling, using 
the Stochastic Series Expansion Quantum Monte Carlo 
method. The main effects we observe are (i) an inhomo- 
geneous distribution of staggered moments, (ii) a reduc- 
tion of the local order parameter with increasing number 
of neighbors, (iii) strong long-range effects, (iv) evidence 
for hierarchical nature of the magnetic ground state, 
and (v) a magnetic superstructure in the spin struc- 
ture factor. These predictions should be tested by high- 
resolution NMR and neutron scattering measurements. 
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